Atomic diffraction from nanostructured optical potentials 
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We develop a versatile theoretical approach to the study of ultracold atom diffractive scattering 
from light-field gratings by combining calculations of the optical near field generated by evanescent 
waves close to the surface of periodic nanostructured arrays together with advanced atom wavepacket 
propagation techniques. Nanometric 1-D and 2-D arrays with subwavelength periodicity deposited 
on a transparent surface and optically coupled to an evanescent wave source exhibit intensity and 
polarisation gradients on the length scale of the object and can produce strong near-field periodic 
modulation in the optical potential above the structure. As a specific and experimentally practical 
example we calculate the diffraction of cold Cs atoms dropped onto a periodic optical potential 
crafted from a 2-D nanostructure array. For an "out of plane" configuration we calculate a wide 
diffraction angle(~ 2 degrees) and about 60% of the initial atom flux in diffraction orders ±1, an 
encouraging result for future experiments. 

PACS numbers: 32.80.Lg, 32.80.Pj, 32.80.Qk, 39.20.+q, 39.25.+k 



I. INTRODUCTION 

The prospect of manipulating neutral atoms and 
molecules by light forces acting at nanometer scale 
lengths offers fascinating but experimentally chal- 
lenging possibilities in many areas of atomic physics. 
Atom optics [0 , atom nanolithography (|, [| , and atom 
interferometry are prime examples. The use of 
high-refractive-index dielectric or metallic nanometric 
objects to produce subwavelength localized light field 
distributions @, f| §J raises the possibility of " inte- 
grated atom optics" in which atoms or molecules can 
be confined, guided or diffracted above nanostruc- 
tured surfaces fabricated to a designed shape |?| ||, [| . 
Tailoring optical potentials for atom control, and 
possibly Bose-Einstein condensates, is analogous to 
the use of micromagnetic fields for a similar pur- 
pose @ 0, H || If 0, 0, £§. Atom diffraction 
from a transmission optical grating provided one of the 
early examples of light-field atom manipulation [fl9| . 
A proposal for |2(| and realization of |2l| a reflection 
light-field grating quickly followed. These early de- 
velopments stimulated many subsequent experimental 
and theoretical studies, and a pertinent review has 
recently appeared p3] . Until now reflection gratings for 
matter waves have been implemented by the formation 



of one dimensional evanescent standing waves produced 
by counterpropagating laser beams undergoing total 
internal reflection through a glass prism. 

The approach we present here contrasts markedly 
with this earlier work on standing waves of sinusoidal 
form. We study diffractive scattering of cold atoms 
from an evanescent field, spatially modulated by an 
array of nanometric objects with high index of re- 
fraction and subwavelength periodicity deposited on 
a glass surface |23|, These evanescent fields and 

their interaction with atoms will exhibit several novel 
features. First, rather than a pure sinusoidal evanescent 
standing wave, nanostructured periodic corrugation 
generates fields containing higher-order harmonics and 
an intricate polarization distribution. Second, spatial 
gradients of field intensity and polarization at length 
scales well below the diffraction limit interact with the 
external and internal atomic degrees of freedom in ways 
that strongly depend on the geometry and material of 
the nanostructures employed. Third, the intensity and 
polarization map above the nanostructure array also 
depends strongly on the intensity and polarization of 
the exciting light source. 

Since the cold-atom de Broglie wavelength is not much 
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FIG. 1: Schematic diagram of cold-atom diffraction off an op- 
tical potential. The atom impinges on the periodic potential 
with incident wave vector ki and diffracts into orders k±i, ko. 
The periodic optical near-field above the nanostructured array 
is generated by light with incident wave vector ko , totally 
internally reflected into ko - The nanostructures are cubes 
of high-refractive-index material with subwavelength dimen- 
sions and periodicity (see text for model details). Note that 
the incident plane of the light Pl is along the diagonal of the 
cubes while the incident plane of the atoms Pa is shifted 45 
degrees and parallel to the cube edges. Note also that the 
diffracted orders k±i are "out-of-plane" with respect to the 
plane of incidence of the light. 



smaller than the characteristic scale length of the optical 
field, an accurate description of atomic motion calls for 
a quantum treatment of external as well as internal 
degrees of freedom. In order to analyze the dynamics of 
cold atoms scattering off periodic optical potentials we 
need to marry two well-developed numerical techniques: 
calculation of the optical near-field and atom wave 
packet propagation. The situation we consider is shown 
m Fig.§ First, we calculate the three-dimensional elec- 
tric field and polarization distribution in the near field 
surrounding the nanostructures. At frequency detunings 
far from resonance, where absorption is negligible, 
this information is used to construct a conservative 
potential for a full three-dimensional treatment of the 
center-of-mass motion of the atom, including ground 
internal states. Using the light field information above 
the nanostructures, we calculate the three-dimensional 
atom-field interaction potential (2^, |2^, |3(J . Since the 
field polarization will vary significantly over the length 
scale of the nanostructures pa ], and since various field 
polarizations may lead to population transfer among 
atomic internal states [ p2[ , we include the internal 
Zeeman states of a 2 S\/2 atomic ground level in the 
calculation of the atom-optical field scattering. Second, 
we apply a time-dependent wave packet method to 
describe the scattering problem of cold atoms diffracting 



from an optical grating with subwavelength periodicity. 
Inclusion of the ground-level internal states leads to a 
three-dimensional, coupled channel problem that we 
solve with wave packet propagation techniques already 
successfully applied to atom or molecule surface scatter- 
ing and quantum molecular dynamics in several degrees 

"0HI H© @- 



33 



35 



of freedom |3fJ, 

After the wave packet representing the atom reflects and 
diffracts from the optical potential, it is projected onto 
final scattering states to yield the desired diffraction 
probabilities j38j. 

In this paper we emphasize general principles and 
seek to establish a methodology without restriction to 
any specific experimental set-up. In fact the numerical 
solution of Maxwell's equations in the near field and 
Schrocdinger's equation for the atomic motion can 
easily be adapted to explore atom manipulation in 
sub- wavelength optical light fields of arbitrary geometry. 
However to illustrate this methodology we present 
calculations obtained with realistic parameters corre- 
sponding to an experiment using a flux of cold 2 5*1/2 
atoms incident at about 40 degrees from the normal and 
scattering off a two-dimensional optical grating with 
subwavelength period. The results of our calculation 
yields out-of-plane diffraction into a few orders, with 
markedly large diffraction angles. The results are 
encouraging for planned experiments. 

The rest of the paper is organized as follows: Section 
II gives details of the calculation of the scattering prob- 
abilities. The sub- wavelength light field calculations are 
reported elsewhere [^8|, so only the essential steps are 
summarized here. In Section III we present a model con- 
figuration and in Section IV we show the numerical re- 
sults for illustration. 



II. THEORY 
A. The optical near field 

The evanescent wave field created by the totally inter- 
nally reflected laser beam is strongly modified by the pe- 
riodic nanostructures on the glass substrate surface. The 
calculation of this intricate electric vector field above the 
nanostructures can be carried out using several methods 
that solve Maxwell's equations in these nontrivial geome- 
tries. These methods are well established in the field 
of scanning near- field optical microscopy (SNOM) [j§ |. 
Among the most widely used are the finite-differences 
scheme |44| , the differential theory of gratings (DTG) 
( H| |||, £7J) and the Green's function method 0. This 
latter approach is well-adapted to study single, finite- 
size nano-objects, and has been used extensively to study 
a wide range of nanostructures p5| Efjj. However, in the 
case of periodic surface structures, it becomes inefficient, 
since it does not explicitly take the periodicity into ac- 
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count. In such cases, the DTG method is more appropri- 
ate, because it solves Maxwell's equations by means of a 
Fourier expansion of each held component. 

To be more specific, we dehne r = (x, y, z) with z taken 
as the direction perpendicular to the surface. For conve- 
nience, we dehne a vector 1 = [x, y) in the plane of the 
substrate surface. We will denote the electric held by 1 
in the glass half-space and by 2 in the vacuum half-space, 
ni being the glass index and the vacuum index. The 
two-dimensional periodicity of the nanostructures defines 
a unit cell of length L x and L y along the x and y direc- 
tions, respectively. The direction of the incoming laser is 
given by {k%,k v L ,k^), with fcj, = n\l^j\ , Ao being the 
laser wavelength in vacuum. 

In the DTG method, an index-modulated zone, la- 
belled 3 and characterized by 713(1"), the high index ma- 
terial, is introduced at the vacuum/glass substrate in- 
terface. We have then to solve in all space Hclmoltz's 



component of the wave vector j„ 
equation: 



obeys the dispersion 
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FIG. 2: Schematic of the three optical regions to which the 
differential method of gratings is applied. 

equation: 

AE W (r) + nf fcgE w (r) = 0; I = 1, 2, 3 

Since the system is periodic along x and y we can ex- 
pand the dielectric constant (i.e. n^) in the modulated 
zone: 



n 3 2 ( r ) = ^2 a mn exp(2l7r( 



nx my 
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(kL) 2 + (7mn) 2 = (klf + (klf + (kl) 2 (3) 



for which we have two solutions, corresponding to "ris- 
ing" and "descending" waves. The rising zero order in 
the glass substrate corresponds to the incident laser held. 
Then, if we know the Fourier components E TOn and B m „ 
of the electric and magnetic field in a plane just above 
the nano-objects, the field distribution anywhere above 
the structures are determined by Eqs.ffl,2. From Eq. g, 
it may be seen that the coefficient 7 mn is either real or 
purely imaginary. The real values of j mn correspond 
to radiative harmonics while imaginary values introduce 
evanescent components. 

The six components of the electromagnetic field 
E (,) (r), B (/) (r) are then deduced from two indepen- 
dent parameters, usually named the principal compo- 
nents, which in the present case are chosen to be the 
y-components E^ y (r) and B® y (r). In order to cal- 
culate the Fourier components of the electric held just 
above the nanostructures, we have to solve a system of 
linear differential equations for E^ y (r) and B"' v (r) in 
the region of space where the index is modulated by the 
nanostructures. This system mixes all the Fourier orders 
of the electric and the magnetic field through the product 
n^E^ y (r) in Helmholtz's equation. Then, for numeri- 
cal applications the Fourier expansion converges at some 
sufficiently large value TV, i.e. the order n and m will 
both vary from —TV to +TV. We solve this system using 
standard bounding conditions for this problem, i.e. we 
assume there is no field arriving from infinity in the vac- 
uum, and the only Fourier components of the incident 
field in the glass substrate are those of zero-order. Fi- 
nally, we can express the Fourier components Emn V and 
Bmn in a horizontal plane just above the nano-objects as 
a linear combination of the y-components of the incident 
fields, Ei,Bi, which define the polarization state of the 
incoming laser : 



where a mn are the Fourier coefficient of nj 2 . It is clear 
that electric and magnetic fields outside the modulated 
zone will also be periodic and thus conveniently written, 



F (2) V _ n-EE F {1)V 



mn i 
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(4) 
(5) 



E©(r) = £ E« e *We*L J (1) 



m.n— — 00 



B«(r) = £ BW e W e *L.l. ; = 1)2 (2) 

m,n— — 00 

In this expression, k|^ n denote the x,y components of 
the wave vector of the field diffracted in the (to, n) order, 

and is given by: k} nn = (k T L + ^,k y L + 3p) . The z 



n-EB n-BE J 

mn J mn J mn ^ 



These transmission coefficients, Tmn 

depend only on the geometry of the sample, the 
frequency and the angle of incidence of the illuminating 
laser. The setup for the field calculation is shown in 
Fig. 2. A more detailed description of this calculation 
can be found in Refs.[^5[ ^?J. With the optical field 
mapping in hand we turn our attention to atom scatter- 
ing in the presence of these fields. 
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B. Calculation of atom-surface scattering 

In the limit of low saturation and a blue detuning 8 
which is large compared to the Doppler shift and natu- 
ral linewidth, the atom-field interaction can be treated 
within the framework of coherent atomic motion where 
spontaneous emission is neglected [^2[ [30|. We will 
consider an atom transition dipole, typical of the first 
alkali 2 S 1 / 2 — > P transition, neglecting fine and hyper- 
fine structure in the excited state, but including the two- 
component angular momentum degeneracy of the ground 
state. The excited level can be eliminated adiabatically, 
which results in atomic motion that is described by a two- 
component, three-dimensional wavepacket that evolves 



on the ground-state manifold of the Zeeman sublevels 
iti 2 ' = T r * m .+ Y, v m 3 m>(r) * m j(r)(6) 



m ' 3 =±2 



In this expression, the operator of the kinetic energy is 
simply given by 



T = V 

2M 1 



(7) 



with M being the mass of the atom. The potential within 
the low saturation limit can be written as p^, [2^, 30 



Vm^ir) = — 22 E q( r ) E q'( r )tig, m 3; l ,<l\je,m e )(j e 



\3g, 



1,5') 



(8) 



In this expression the terms in parenthesis are Clebsch- 
Gordan coefficients with j g and j e being the total an- 
gular momentum of the ground and excited state re- 
spectively. The reduced dipole moment is denoted by 
d and the field enters through its spherical components 
E q (r), q = 0, ±1. This expression forms the basis 
of numerous studies of atomic diffraction from stand- 
ing evanescent waves as has been used by many authors 
" |2T 



H |22|, |29j, |30|, 43 . In general the polarization state of 
the electric field above the nanostructured objects is a 
very complicated function of space. Hence, according to 
Eqs.||,||, transitions between the different Zeeman sub- 
levels can occur during the interaction of the atom with 
the light fields. As a consequence, we performed the col- 
lision on the coupled ground state levels. 

One way of proceeding is to expand the wavefunction 
into plane waves parallel to the surface which results in a 
set of coupled diffraction channels for the direction per- 
pendicular to the surface . Diffraction then occurs as 
transitions between these diffractive states, the probabil- 
ities of which can be calculated by a semi-classical treat- 
ment within the Landau-Zener theory, p2]| or by numer- 
ical wave packet propagation fi~3| . This treatment cor- 
responds to the well-known CCWP (close-coupling wave 
packet) method used in standard atom or molecule sur- 
face scattering |42| . However, with increasing diffraction 
orders one may need to take a large number of diffrac- 
tive states into account, together with all the coupling 
matrix elements among them. In fact it is only recently 
that modern and very efficient numerical wave packet 
propagation techniques have permitted the calculation 
of diffractive scattering by solving Eq. 6 directly on grids 
in real space [M], |3^, |{J . Briefly, one simulates the col- 
lisional process by a wave packet propagation from the 
initial state to the final interaction free zone, where it is 



projected onto any desired observable, which in our case 
is the total population in the different diffraction chan- 
nels. To achieve a good energy resolution, one needs to 
construct an initial wave packet that is sufficiently large 
such that its energy width is negligible. From a numer- 
ical point of view, this is very disadvantageous, since it 
requires a large grid in the direction perpendicular to 
the surface. Therefore, we used a method proposed by 
Mowrey and Kouri [02] , who started with a spatially nar- 
row wave packet that comprised a wide range of energies, 
and extracted energetically resolved results by projecting 
onto asymptotic states of well-defined energy, E . 

The initial state of the atom is taken to to be a specific 
Zeeman sublevel m° of the ground state manifold with 
center of mass motion described by the box-normalized 
wavefunction 



^o(r) = (2neL 2 x Ll) 3 x 

exp (-0 - z ) 2 /4£ 2 + ikfz + ik, • 1) 



(9) 



with the initial transverse momentum = (kf,kf) and 
the initial momentum in z-direction described by a Gaus- 
sian distribution centered around kf. The propagation 
of this three dimensional wavepacket is performed on the 
coupled surfaces defined by Eqs. until the final scat- 
tered wave function . (r) is entirely in the asymptotic 
region. 

The (unnormalized) modulus square of the transition 
amplitudes is given by projection of the final wave packet 
onto diffractive states WA : 



IC"(£)| 2 = 



dr 



-(k,+G„ 



(10) 
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FIG. 3: Top panel: schematic view of the 100 nm TiCh cubes 
deposited on a silica surface. Bottom panel: Light field inten- 
sity distribution 125 nm above the nanostructures (the arrow 
kJJ indicates the plane of the incident laser of mode TM and 
k[' the plane of incidence of the atoms. The contour lines 
correspond to (in units of |-E|;„ C ): 0.397, 0.412 (dashed lines) 
and 0.427, 0.442 (full lines) 
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) denotes the reciprocal lattice 



where G mn 

vector and h z mn is determined by energy conservation: 



¥ mn = y/2M£ - (ki + G r 



(11) 



The probability to diffract into a specific order (m,n) as a 
function of total collision energy is then given by M, H2J : 



Pmn (£) 



— |_i inn 



This quantity will be calculated for a selected set of 
parameters and discussed in the next section. For 
the three-dimensional, two-component wave packet 
propagation we used the FFT-Split-Operator scheme, 
which in the field of quantum molecular dynamics or 
molecule-surface collisions has proven to be a fast, stable 
and efficient method to perform this task p8[ . Briefly, 
in this method the total quantum mechanical short-time 
propagator is approximated by 



exp 



.At 



T r + Vi 
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V i i 
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(r) 
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T r + V_x _i(r) 



(13) 
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i(r) 



exp 



At { T r 



2^ V T r 



This form of symmetric splitting is correct through 
second order; it is unitary by construction and thus en- 
sures numerical stability |49j. The operator acts on 
the two-component wavefunction to perform a short-time 
propagation from time t to t + At. The action of the ki- 
netic operator is calculated in Fourier space, where it is 
a simple multiplication with a phase factor, and so ev- 
ery time step the two components of the wavefunctions 



need to be Fourier transformed which can be done very 
efficiently with three-dimensional FFT algorithms. With 
this technique, the initial wavepacket Eq.(^|) is propa- 
gated until it is entirely in the asymptotic region. From 
which distance the wavepacket can be considered to be ef- 
fectively free has to be checked carefully since it depends 
on many parameters. This is especially important in the 
field of cold collisions where low energies are considered. 
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Even though not applied in the current approach, we 
note that the numerical effort can still be reduced by an 
adiabatic correction of the initial state |3j| |fj|, [50| , an- 
alyzing the flux out of the scattering region instead of 
projecting onto final diffractive states |30|, |l) or using a 
filtcr-diagonalization scheme |p2| . 



III. NUMERICAL IMPLEMENTATION 

A. Model Configuration 

We model a typical experimental setup with a regular 
square lattice of TiC>2 cubes (100 nm on a side, index 
of refraction n — 2.1) deposited onto a flat silica surface 
(index of refraction n = 1.5) at a center-to-center dis- 
tance of 250 nm (Figs.0J|^). These structures are illu- 
minated by an evanescent light field created by an incom- 
ing laser with a vacuum wavelength of 850 nm, intensity 
of 80 W/cm 2 , and subject to total internal reflection. 

The plane of incidence of the laser beam (denoted kg 
in Figs. 0J||)is chosen to be diagonal with respect to the 
rectangular nanostructured pattern, and its polarization 
is taken to be within the plane of incidence (TM polariza- 
tion). The angle of incidence with respect to the surface 
normal was chosen to be 60 degrees. Under these condi- 
tions, the flat surface without the nanostructures would 
give rise to an evanescent wave with a decay length 700 of 
about 170 nm, and the atoms would not approach closer 
than about 200 nm to the surface. Hence the influence 
of the attractive van der Waals atom-surface potential 
can be neglected. For different parameters however, it 
might become important |^7j , and can easily be included 
in Eqs.|6|H. The Fourier decomposition Eq. [I] required 
N = 10 terms for convergence. 



B. Diffractive scattering 

We have chosen a model atom with ground state 2 Si/2, 
the mass of atomic Cs, and transition dipole moment cor- 
responding to the 6 2 Si/2 — * 6 2 P 3 / 2 atomic Cs transition. 
The initial internal state was taken to be rrtj = —1/2 
and the initial Gaussian distribution of the perpendicu- 
lar motion was centered around k\ = 0.87 nm -1 with a 
width parameter £ = 5.00 nm. The inital transverse mo- 
mentum is taken to be ki = 0.73 nm -1 , ki = 0.0. This 
corresponds approximately to cold Cs atoms produced in 
a MOT after a free fall of about 1.5 cm before colliding 
with the nanostructured surface inclined at an angle of 
about 40 degrees between the vertical axis and the sur- 
face normal. Under these conditions, the incoming and 
specular direction are in the y — plane, as indicated in 
Fig. 0. We found converged results for a cut-off distance 
of about 600 nm. The asymptotic wavefunction is then 
analysed to yield the diffraction probabilities P mn (£) as 
outlined in the previous section. In all calculations, grids 
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FIG. 4: Light field intensity distribution in the y, z plane 
together with a schematic view of the Ti02 nanostructures 
deposited on the silica surface. The contour lines correspond 
to (in units of \E\ 2 inc ): 0.080, 0.135, 0.189, 0.241, 0.293, 0.346, 
0.400 (dashed lines, from top to bottom). 



of 16 points in the x and 64 in the y direction were used 
and 1024 points in the ^-direction. 



IV. RESULTS 

A. Sub-wavelength optical near-field 

In Fig.|| we show the geometry of the Ti02 cubes (a) 
together with the intensity distribution (b) at a distance 
of 125 nm above the surface. The arrow indicating the 
plane of incidence of the laser beam is denoted by kjj and 
plane of incidence of the atoms by k^. One can clearly 
see that the light fields bear the periodicity of the nanos- 
tructures. We found that illuminating the nanostruc- 
tures with the plane of laser incidence aligned along the 
cube diagonals rather than along the sides yields peri- 
odic potentials with steeper gradients (more pronounced 
localization) even at distances far above the surface. At 
the illustrated distance of 125 nm above the surface, 
we still find pronounced periodic field intensity modu- 
lation. This strong localization of light intensity above 
high-refractive-index structures is a well-known charac- 
teristic of the TM illumination mode. It has been the- 
oretically modelled and experimentally verified by the 
SNOM technique for a number of nanostructured ob- 
jects j|. In Fig. |] we show a cross section of the field 
intensity distribution in the y, z plane above the nanos- 
tructures up to 400 nm. One clearly sees how the inten- 
sity contrast above the nanostructures decays with in- 
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FIG. 5: Distribution of the spherical components of the elec- 
tric field 125 nm above the top of the nanostructures. The 
contour lines correspond to (in units of |P|;„, C ): (a) 0.377, 
0.387 (dashed lines) and 0.397, 0.407 (full lines): (b) 0.020, 
0.025 (dashed) and 0.030, 0.035 (full lines): (c) 2.3- 10" 4 , 
4.6-10" 4 (dashed) and 6.9-10" 4 , 9.2-10" 4 (full lines). 



creasing distance. Depending on the kinetic energy, the 
falling atoms will penetrate to different heights above 
the surface, therefore experiencing a different amplitude 
modulation (contrast) at the plane of the mean classi- 
cal turning point. It is important to note that even if 
the contrast diminishes with distance, the periodicity re- 
mains the same; and that even if the atoms are diffracted 
at large distances from the surface, the diffraction angle 
will still be controlled by the ratio of the atom de Broglie 
wavelength to the optical grating period. In order to rig- 
orously establish the optical potential that governs the 
atomic motion, we need the full information of the field 
above the surface, i.e. the three spherical components of 
the electric field E+ , E_ and Eq , as can be seen from Eq. 
|8|. The quantization axis was chosen to be perpendicu- 
lar to the plane of laser incidence. From Fig. [3| one can 
see that the different components of the light fields are a 
complicated function of space that will interact with the 
multi-level internal structure of the colliding atom. Note 
that the Eq component is one order of magnitude smaller 
than E + or P_ . 



B. Diffraction probabilities including ground state 
degeneracy 



In this section, we show the results for diffractive scat- 
tering, using the light field parameters as detailed in 
the previous section. We stress that these probabilities 
are calculated using the Cs 6 2 S , 1 / 2 ~ * 6 2 P 3 / 2 transition 
dipole and atomic mass, but that no other specific atomic 
parameters were used. This particle is still a model atom, 
however, since we have ignored hyperfine structure in the 
ground and excited states. The center-of-mass motion is 
treated entirely quantum mechanically, and the polar- 
ization state of the electric field with its spatial varia- 
tion rigorously included in the calculations. Since we are 
considering a 2-D potential surface, the diffraction takes 
place in two spatial directions labelled G mn , which are 
the reciprocal lattice vectors corresponding to the x— and 
y— directions, respectively. The probability P(m, n) to 
diffract into a given order (m, n) as a function of energy 
is given by Eq. 12. Figure plots these diffraction prob- 
abilities P(m, 0) (top panel) and P(0,n) (bottom panel) 
as a function of total energy. The symmetrical results for 
negative values of m,n are not shown for clarity. With 
the definitions employed (m, 0) corresponds to diffraction 
within the plane of incidence of the atoms while (0, n) 
corresponds to diffraction perpendicular to this plane. 
These two cases will be called 'in- plane' and 'out-of- 
plane' diffraction, respectively. Comparing the top and 
bottom panels of Fig.||, we see that probabilities for in- 
plane and out-of-plane diffraction differ dramatically. 

The top panel of Fig.^ shows the in-plane diffraction 
probability as a function of energy for fixed initial trans- 
verse momentum of kf = 0.73 nm -1 . For low energies, 
we find negligible diffraction. Most of the final popu- 
lation is in the (0,0) channel which represents specular 
reflection. As the energy E increases, P(0, 0) decreases 
and falls off to almost zero for energies greater than 2.5 
mK (E = /cbT). The higher orders of in-plane diffraction 
P(m, 0), to = ±1,±2,±3,.. are not significantly popu- 
lated either. The bottom panel of Fig.^ shows where 
the scattering flux has gone. Almost the entire popula- 
tion is found in the out-of-plane diffraction orders (0,n), 
77 = ±1,±2, ±3. The fact that we calculate a rapidly 
decreasing in-plane diffraction corresponds to the well- 
known averaging of the transverse scattering amplitude 
of the atomic motion parallel the diffracting structures 
p2j. With our parameters, we find an interaction time 
of about 3 /zsec, which together with the initial trans- 
verse momentum of kf = 0.73 nm" 1 implies that the 
atoms "sample" the modulated potential over a trans- 
verse distance of about 800 nm, greater than the peri- 
odicity of the potential by more than a factor of three. 
For the out-of-plane scattering there is no such averag- 
ing and the situation is analogous to in-plane diffraction 
at normal incidence. It is worthwhile noting that this 
out-of-plane diffraction at grazing incidence (in a quite 
different experimental arrangement than the one envis- 
aged here) has been observed experimentally 153] . 
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collision energy (mK) 



FIG. 6: Upper panel: In-plane diffraction probability P(m, 0) 
(the plane containing ki), Lower panel: Out-of-plane diffrac- 
tion probability P(0, n) (perpendicular to the plane contain- 
ing ki), as a function of total collision energy(_E = fcsT). The 
arrows indicate the cases discussed in the text. 

Using the results shown in Fig.^j, one can now sim- 
ply read off the probability for specific collision ener- 
gies. At an energy of 2.35 mK, which corresponds to 
an angle of incidence of about 40 degrees (indicated by 
the left arrows in the top and bottom panels of Fig.^), 
we have almost the entire population equally distributed 
in the three diffraction orders (0,0) and (±1,0). Fig- 
ure [7] shows the complete set of diffraction probabilities 
P(m, n) for this case. One sees clearly that the only sig- 
nificant diffraction is out-of-plane, and the three domi- 
nating diffraction channels have almost equal probability 
of 30 percent each. The diffraction angle, given roughly 
by the ratio of deBroglie wavelength to grating periodic- 
ity, is quite large-about 2.6 degrees. This geometry could 
thus be a very attractive for the realization of an atomic 
beam splitter or for atomic interferometry. 

As a second example, we have shown in Fig.|| the 
situation at slightly higher collision energy of 2.52 
mK. With the chosen initial transverse momentum 
this corresponds to an angle of incidence of about 38 
degrees. The maximum of the final state population is 
now concentrated in the orders (0,±1) while the specular 
direction is suppressed to about 10 percent. At this 
slightly higher energy, we find also the diffraction peaks 
(0,±2) to be populated. Still, the ratio between P(0, 0) 
(or P(0, ±2)) and P(0, ±1) is about 3, and the angle 



FIG. 7: Diffraction probabilities P(m,n) for Cs atoms with a 
total energy of 2.35 mK corresponding to a drop height of 15 
mm (left arrows in Fig. 3) and an angle of incidence (between 
the vertical axis and the surface normal) of 40 degrees. The 
diffraction angle between the orders (0,±1) is ~ 2.6 degrees 



spanned by the two main diffraction peaks P(0, 1) and 
P(0, —1) is about 2 degrees. 



C. Diffraction probabilities neglecting ground state 
degeneracy 

In the previous section, we have calculated the diffrac- 
tion probability including the ground state degeneracy. 
As can be seen from Eq.|8|, the intricate spatial distribu- 
tions of polarization in general do not allow treatment 
of the problem as a one-surface scattering event, but re- 
quire the solution of multiple internal states coupled to 
the optical potential. However, as one sees from Fig. 3, 
the electric field is dominated by the E + and _E_ compo- 
nents, with the Eq one order of magnitude smaller (the 
axis for the polarization state is chosen to be orthogonal 
to the plane of incidence of the laser beam). One sees 
from Eq.|| that the coupling between the ground state 
sublevels in the present case of an S%/2 ~ * P3/2 transi- 
tion is due to the Eq component. Hence if the component 
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tational burden, can yield satisfactory results. Thus in 
this particular case a one-surface calculation gives very 
good agreement with the full results, even for the weakly 
populated diffraction orders. In general, however, with 
atoms having higher ground-state multiplicity and with 
different optical materials and excitation geometries, a 
full coupled-surface calculation will be necessary. 




P(m,n) 
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FIG. 8: Same as Fig.fj], but for a total energy of 2.56 mK, 
corresponding to a drop height of 16 mm, (right arrows in 
Fig.^|) and an angle of incidence of 38 degrees. The diffraction 
angle between the orders (0, ±)is ~ 2.0 degrees. 



V. CONCLUSIONS 

In this paper we have investigated the diffraction of 
cold atoms by highly structured subwavelength optical 
potentials generated from evanescent fields. Our ap- 
proach includes a three dimensional quantum treatment 
of the atomic center of mass motion. We take into ac- 
count the spatial distribution of rapidly varying polar- 
ization states of the nanostructured optical fields and in- 
clude the effect of these polarization changes on the atom 
ground internal state populations. 

As an illustration, we have chosen a model system that 
corresponds to cold cesium atoms (without nuclear spin) 
diffracting from a nanostructured surface illuminated un- 
der conditions of total internal reflection. The interaction 
of cold atoms with these light fields is calculated in the 
limit of large detuning and negligible absorption. For ex- 
perimentally realistic initial conditions, we find diffrac- 
tion angles on the order of 2 degrees with about two 
thirds of the initial atomic flux concentrated in the first 
two diffraction orders. These structures may therefore 
prove useful in wide-angle atomic interferometers. 
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